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Abstract 

We consider the heat-kernel on a manifold whose boundary is piecewise smooth. 
The set of independent geometrical quantities required to construct an expression 
for the contribution of the boundary discontinuities to the Ci heat-kernel coefficient 
is derived in the case of a scalar field with Dirichlet and Robin boundary condi- 
tions. The coefficient is then determined using conformal symmetry and evaluation 
on some specific manifolds. For the Robin case a perturbation technique is also 
developed and employed. The contributions to the smeared heat-kernel coefficient 
and cocycle function are calculated. Some incomplete results for spinor fields with 
mixed conditions are also presented. 
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1 Introduction 



In the theory of quantum fields on a curved background space, the short-time ex- 
pansion of the heat-kernel has been found to be of paramount importance. The 
coefficients in the expansion are involved in, for example, statistical field theory and 
the calculation of vacuum energies. 

Much effort has been expended in finding explicit expressions for the coefficients 
but, in view of the rapidly increasing complexity, the evaluation of higher and higher 
coefficients would seem to reach a point of diminishing returns. Rather than pursue 
this path, the present work seeks to consider a situation where the domain itself 
is generalized. We do this within the context of the C2 coefficient which is of 
particular importance in 4 dimensions being linked to the 1-loop quantum effect on 
the Einstein field equations, as well as to the conformal anomaly, of relevance to 
quantum cosmology. 

So far this coefficient has been determined in the case where the manifold is 
closed U or has a smooth boundary In this paper, we determine C 2 when 

the boundary is piecewise smooth. This is an extension of our previous work on C3/2 
in the piecewise smooth case Q. Our approach is to find the most general possible 
expression for the contribution to C2 arising from the discontinuity - a non-trivial 
problem in itself - and then to confine it using the required conformal symmetries 
and special-case evaluation. 

This knowledge then enables us to find the expression for the change in the 
effective action under a conformal transformation. The latter provides a follow-up 
to previous calculations of the effective action on specific manifolds with smooth 
boundaries § §, § % §. 

The physical motivation for this analysis is that piecewise smooth boundaries, 
and manifolds, occur in various idealized situations; for example in simplicial ap- 
proximations to general relativity and to quantum gravity, || |TU|] . 

Also, internal vertices (conical singularities), which are the periodic versions 
of the present structure, appear in the theory of quantum fields on a black hole 
background and have been the subject of some activity. 
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2 Background 

The quantity of interest is the term in the small-t expansion of the integrated 
heat kernel on a d-manifold At with positive-definite metric 

m = j M avsw^w = ^J v *' ID 

for the massless scalar field operator 

A = -V 2 + ^R (2) 

on At, where the fields are real, and obey Dirichlet conditions at the boundary 
dAi. £ is a variable coefficient of coupling to the Ricci scalar R. A is conformally 
covariant if £ is equal to 

d-2 



4(d-l) 

Although we consider a massless field for simplicity, our results are easily generalised 
to the non-zero mass case. 

In principle, in singular situations, there is the possibility of \ogt terms in the 
expansion. There are none such in the present situation and they are henceforth 
ignored. 

The heat-kernel coefficients can be written as integrals of local geometrical quan- 
tities such as the curvature tensor over Ai and its boundary dAi if one exists. A 
traditional method of determining the integrand is to write down the most general 
possible expression and then use known values of on particular manifolds to 
confine it. This shall be our approach. 

A convenient method of calculating the heat-kernel coefficients in specific cases 
where the eigenvalues Aj on At are known is via its connection with the generalized 
zeta function 

Cm(*) = EV (3) 

i 

where an analytic continuation is involved for TZe(s) < d/2. Cm{ s ) can be written 
as the Mellin transform of the integrated heat-kernel, yielding its connection with 
theCf JTT|: 



(4) 
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for k = 0, 1/2, 1 • • •. In this way, using the Aj on a particular manifold to calculate 
the zeta function gives us the values of the heat-kernel coefficients on that manifold, 
and hence information about the general expressions for the coefficients. 

If A and the boundary conditions are conformally covariant, then from the 
change in the Aj brought about by a conformal transformation, it is easily shown 
that Cm(0)j an d hence C^fX, is conformally invariant. In other words, is con- 
formally invariant in 2k dimensions - this places a large restriction on its form. 
More generally and more quantitatively, the conformal variation of the heat-kernel 
coefficients under a transformation g^ u — > e 25tx) ^g^ y for small Scu is |12 



5Cl d) = (d- 2k)Cl d) [5co(x)} + 2Cl d } 1 [J5co(x)} (5) 

where 

J = (d-l)(£-&) V 2 
and the smeared coefficients C^[f(x)] are defined as 

/ J d x^(x\e- tA \x)f(x) = — ^ £ 3f°t/(*)]** ( 6 ) 



fc=0,l/2, 



Although our main interest is in determining the effect of boundary discontinu- 
ities on C?2 , we give the result for the smooth boundary case first. It is found that 



C {d) = — I d d xJ~g 
2 360 Jm V 



— / d^xVh 
560 JdM 



2R^R pvpu - IBTR^ + 5(6£ - l) 2 i? 2 

320 , 88 2 40 o „„„„ 

tm 3 ntiK 2 H k — AR^Kw 

21 7 21 



360 

-AnR^n^nv + l6R^ pa n^n p n ua + 10(1 - 6f) (2Rk - 3V n R) 



(7) 



is the inward-pointing normal to the boundary. With the induced metric 
h^v = (jfj,v — n^riy we then define the extrinsic curvature tensor = K vtl = 
—h^h^V a n/3, with k = k^, tr/t 2 = k^k^ . In our conventions k is positive on 
the surface of a ball. As we would expect from (ED, (|?p is conformally invariant in 4 
dimensions if £ = 1/6. 

We note that the coefficients in (0) are independent of the dimension. This is a 
general feature, and can be proved by multiplying Ai by a circle. We refer to the 
book by Gilkey |H| for general information regarding heat-kernel asymptotics. 
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3 Boundary discontinuities 



We shall now consider the situation where DM. is not smooth, but is made up of a 
number of boundary parts, dM.{. The intersection X^ between two adjacent such 
parts is a manifold of dimension d — 2 (i.e. codimension 2). We also require that 
the Xij be closed and smooth. The hemiball and ball x interval are examples of 
this situation. In general the spectral coefficients will now contain "edge" terms 
involving integrals over the X^, as well as volume and boundary terms. 

In addition to the boundary normals nf , we define normals nf^^ to 2^ pointing 
into dAii. Xij has induced metric 7^P- Then on X^, suppressing the (ij) index for 
convenience: 

lnu — — n iX n v = h ] ^ v — n^nl 

h U = 9^ - ri^nl (8) 
For a dihedral angle 0y between dM.i and dA4j, the normals have the relationship 

rij = nf cos 9 + nf sin 9 

rij = nf sin 9 — nf cos 9 (9) 

We limit ourselves to the case where 9 is a constant on each intersection. 
The codimension-2 contribution to will have the form 

W d d ~ 2 x^S i3 

(ij) ^ 

for some scalar depending on the local geometry. Our first restriction on is 
that Sij — > a~ 2 Sij as g^ v — > a 2 g^ u , where a is a constant, to make cf 1 invariant 
under a global conformal (scale) transformation. Expressions with this property can 
be split into two types: intrinsic and extrinsic. Intrinsic terms must be first order 
in the Riemann tensor, contracted with the various metrics and normals. Extrinsic 
terms do not involve the Riemann tensor, but contain two V M operators which act 
on the normals and metrics to produce extrinsic curvatures. It is easily shown using 
dimensional arguments that these are the only types of term possible. 

Intrinsic terms: 
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With intrinsic terms, the objects we have to work with are, on each intersection: 

Rpvpa Rpupa Rjivpv Rpvpa 

9pu n.* n£ raj, 

R\ V pa an d Rpvpcr are the Riemann curvature tensors of dM.i and Ijj respectively, 
formed using the appropriate induced metrics together with the projected covariant 
derivative. h l , b? v and 7^ are redundant since from the relations @ they can be 
expressed in terms of other quantities. 

To construct a general form for <Sy, we have to consider, as far as intrinsic terms 
are concerned, every possible independent contraction of the expressions above which 
involves a single Riemann tensor. From @, we never need to mix up different 
indices i and j in any particular term, since any expression which involves a mixture 
of indices can be written as a linear combination of terms which are pure in i or j. 
The possible contractions are 

pi = R p\ = Ri p3 = R 

p\ = Rpurfn" pi = R^nffi pi = Rp U nfn^ 

p\ = R'pvKn" p\ = R^n?n v i n p i nl 

and obviously also p J 2 , pi etc. All others are expressible in terms of the above set 
using the symmetries of the Riemann tensors, or vanish since 

3« 5 ^ 5 CjL 1 ri 

rL pupa ll i lx pvpcr l <'i 1 >'pupcr n i U 

In fact, not all of the above quantities are independent, since Gauss' equations 

DM _ M i _ M K i I ha vfl d/i (in\ 

and (defining kL, to be the extrinsic curvature of lij with respect to dM.,j) 

R^upa = K ip K vo ~ K ia K l>p + Ivlp URiaPi 

provide relations between them, assuming we include the required extrinsic terms, 
which will be dealt with separately. Contracting (|T0| ) and ([TT| ) with gPg ucr , and ( |10D 
with g^n u n eT , we get, respectively 

2pe = Pi - A + ^ - tr^ 
2p\ = p i 2 -p 3 + R 2 i - tin 2 , 
Pi = P\ ~ P7 + fi^ - rf p K i lxv n p i n v i 
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A further relation can be found using (|^): 

p\ = /^cosec 2 ^ - pgCot 2 # + 2p4Cot^ 

We choose the p\ to p 4 as our basis. must be symmetric in % and j, so it has 
the form 

Sij = diR + a 2 \ Ri + Rj^j + a 3 R + a^R^ v (n^n" + n^n v ^j + extrinsic terms (12) 

Note that technically one of the R^n^n^ is still superfluous - we could more eco- 
nomically replace the last term by a^R pv n^n^ . However we choose not to do this 
since the expression we have vanishes on an Einstein space, which will turn out to 
be convenient. 

Extrinsic terms: 

These terms must contain either two extrinsic curvature tensors, or the derivative 
of one. In the former case, we have the combinations 

rv tii/ rv pcT ^pv^po ^pu^pa 
™ pumper ^pu^pa ^pu^pa 

plus the expressions created by swapping i for j. The extrinsic curvature of with 
respect to dA4i is defined as k 1 = 7c* = — 7^7^ V Q n^. We contract the above with 

7 M ^ nf rij nf fij 

It is not necessary to include terms where all the indices of k 1 v or r& are contracted 
with 7, since these can be rewritten in terms of R l and kP : 

7 /i V /3 < 3 = -7 MQ 7" /3 V a n^ = «f cosecfl - Aefcotfl 

where we have used (^). Additionally, ((^) means that we can choose not to contract 
with an n M or n M of different indices i and j. 
We are left with the contractions 



k\ 




Li _ 
hv 2 — 




h 




fc 4 = 




kl 




Li _ 
fi, 6 — 


^pv^pj^i^il 9 " 


7 


= ^p^ v l pa 


Li 

fi, 8 — 






= K^Ur u Y a 


fcio = 
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as well as k{, etc. Everything else vanishes, since 

K pu n i = K /j,u n i = ^]j,u n i = 

The fact that we have restricted 9 to be constant provides an extra constraint 
'y'* L 'V u (nfn J a ) = 0, giving 

<7 W ^ = ^k 4 = -k\ = -k{ 

and we discard /c 4 . 

In another possible representation, we could use only k % and k? , since from (g) 

"T'-f^Kp = -l m l v ^ otfp = 7^7^ (<3cosec0 + w^cotfl) 

However, we prefer some less alien expressions, and will use in our calculations 

k% = k 7 cosec 2 9 - 2k 9 cot9cosec9 + k l 7 cot 2 9 + 2 (fcgcosecfl - fcgcotfl) + k\ 

tin 2 = k J 8 cosec 2 9 - 2k 10 cot9cosec9 + k l 8 cot 2 9 + k\ + 2k l 2 

KiKj = — (k 7 + kjj cosec^cot^ + kg (cosec 2 ^ + cot 2 6^ + (k\ + k^j cosec^ 

- (hi + fcg) cot9 + k 3 
KiRi = k g cosec9 — k l 7 cot9 + k\ 
KiRj = k 7 cosec9 — k$cot9 + k\ 
% i = k 7 
tin 2 = k l 8 

tVKiKj = k±o 

These are manifestly linearly independent (assuming the fc's are). In fact, since 
k\ = k\ we could discard, for example, txK 2 or tin 2 but to keep i «-> j symmetry we 
do not do this. 

Derivative terms: 

It remains to consider terms which contain the derivative of the extrinsic curva- 
ture. The tensors involved are 

V \ K% vp ^V K lp 
^fJ. K up ^H K ip 
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where V M and V M are the covariant derivatives projected onto dM.% and X^- respec- 
tively: 

Any other projections of the covariant derivative give nothing new, e.g. 

n v Xl - k' 1 K iu 

which we have already taken into account. This leaves the possible contractions 
plus 5{ etc. The Codazzi equation 

v^-vx^RppMX (is) 

gives a link between the first two: 
and we discard 5\. 

Since X^- is closed and smooth, any term which can be written as the divergence 
of a vector field V^V^ on X^ can be ignored as its integral over vanishes. The 
only vector field on X^ with the appropriate scaling we can construct is 

Vt = 7"XX 

Then 

^ ' pVi = K i K% p,i/ n i n i ~~ K% pu K i + $1 — $3 

The first two expressions can be written in terms of the fc's, so we can discard S^. 

This concludes our enumeration of the terms involved in the codimension-2 con- 
tribution to C^. We have finally 

Sij = diR + a 2 [Ri + Rj^j + a 3 R + a^R^ (nfn^ + n^n 1 ^ 
+6i (k? + nf) + b 2 (tr^ 2 + tr/cj) + 6 3 (X 2 + £ 2 ) 

+64 (tr£ 2 + tmfj + 65 (/Cj/Cj + KjKi) + Ci (/CiKj + KjKj) 

+C 2 KiKj + C 3 /%£j + C^TKiKj + (Zi (V~Kj + V~./Cj) (14) 

The a, b, c and d coefficients are functions of 0, but not of the dimension. The a's 
multiply intrinsic terms while the fe's and c's, as we shall see, form two naturally 
separate groups when = tt/2. The problem is now to evaluate these 14 coefficients 
by obtaining 14 independent constraints. 
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4 Conformal invariance in 4 dimensions 



A great deal of information about the coefficients can be derived by using conformal 
invariance. Under a conformal transformation — > e 2u} ^g pu , we have 

R„v -> Rnv - g^V 2 uj + (d - 2) (VpuVvU - ^V p cjV p cj - V^V^) (15) 

R% - ~ + (d-3) (V^V^ - ti^VpuNOu - V,V u uj) (16) 

R pv -> ^ - 7m^V 2 cj + (d - 4) (y^fJVuU - 7^ V p uV p uj - V M V^) (17) 

e w «, - ^V ni a;) (18) 

^ - 7/-/ V- J) (19) 



'Si/ 



'V^ — ¥ 

We use the above to first order in uj to work out the change in «Sjj under a 
small conformal transformation g pv — > e 26u} ^g pu . With general £, there is an added 
complication in that is not in fact conformally invariant, and we must use 
equation @. This involves the codimension-2 contribution to [/], which is 
known to be [[5] 

(20) 

The change in C% then turns out to have the form 

(50^ - 2C[ d) [J6u]) I = 0{d - 4) + Y, $ j d d ~ 2 x^ Sij 
-Vd / d d ~ 2 x^f (dfWpSuy, - KiV^.8w + rf u n^5uj u - /%V~6o;) 

f..\ J 1i A 



E/ d d - 2 X ^U^-^\ J8u (21) 



(y) ' 

where J is given beneath equation ([5]). The extra term, which vanishes in the smooth 
boundary case, is equal to the boundary divergence which comes from the conformal 
variation of the volume and boundary parts of C% ■ We keep the coefficient rjd 
general since we wish briefly to discuss spin-1/2 later. 

The right-hand side of (|2T|), which we must set to zero in 4 dimensions, is ex- 
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pressible in terms of the 7 linearly independent quantities 

(«iV ni + KjV nj )5uJ {KiVnj + KjV ni )80J 
(£iV ni + KjS/n^U {KiV nj + KjV ni )SuJ 

«^ + <^;)<^ 

where indices on 8uj denote covariant derivatives on Ai. Any other quantity which 
occurs can be written in terms of these, using (^) and the relations between Lapla- 
cians on Ai, dAii, dAij and Ty. V 2 5u can be ignored since it is a total divergence 
on X. Setting the coefficients of these to zero at d = 4 gives us 7 constraints on the 
a, b, c and d. We find 

(4a 2 + 6&i + 2b 2 ) sin 9 + 2b 5 + (2c x + di) cos = rj 4 cos (22) 

2& 5 cos 9 + 3c 2 sin 9 + 2ci + c?i = r/ 4 (23) 

(4& 3 + 2& 4 ) cos + 3ca sin + 2c 3 + c 4 = (24) 

-6ai - 8a 2 + 46 3 + 2& 4 + 36 5 sin 9 + (2c 3 + c 4 ) cos 9 = (25) 

6ai + 4a 2 (l + cos 2 6) + (2a 4 + 3rf 4 ) sin cos = (26) 

(6ai + 8a 2 ) cos + (2a 4 + 3di) sin 9 = (27) 

3di = r? 4 (28) 

where a± is a correction for non-conformal coupling: 

*- + 5(<-i)P^) ' 29 > 

Also, from conformally varying (0) we find that 

rf " 6 1 fort 

^=-90-=>% = -45 (30) 

5 The smeared coefficient and conformal invari- 
ance in 6 dimensions 

There is a further conformal invariance we can use, in that || 

sc (2k + 2) [f] = Q &s ^ ^ f _^ e _ 2aaJ/ (31) 
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as long as A is conformally covariant in 2k + 2 dimensions. 

To implement this, it is necessary to calculate the smeared coefficient [/] 
using (|5]). We do not give it here since it is rather involved. Applying the conformal 
variation above in 6 dimensions with £ = 1/5 then yields the following independent 
terms in the integrand of the codimension-2 part of 5C^ [/] : 
/(«iV ni + KjV nj )5uj f{K>iV nj + KjV ni )5u 

/(£jV ni + KjV nj )5uJ f{KiVn 3 + KjV ni )8uj 

/(kJX + fV 2 5uo 

(V n JV ni + V nj fV nj )8oj (V n JV nj + V % ./V n JcW 
All other terms can be written in terms of these, up to divergences on Xij. Somewhat 
surprisingly, it turns out that the volume and boundary parts of 5C^[f] contain 
no divergences, so that the codimension-2 part vanishes by itself. We therefore set 
the coefficients of all the above terms to zero. The first seven simply give us the 
equations we have already derived from 4-dimensional conformal invariance. From 
the remaining three, we get, respectively 

6S1 + 8a 2 + 3a 3 = (32) 
-12ai - 16a 2 - 5(&i + b 2 ) sin 2 9 + 863 (l + cos 2 0) 

+46 5 sin 9 + 8c 3 cos 9 + ( 4ci - —d x - — ] cos 9 sin 9 = (33) 

V 2 54/ 

(-12a! - 16a 2 + I663 + 96 5 sin^) cos 9 
+4c 3 + (9c x + 5c 2 sin 9 - — - 3d^ sin 9 = (34) 

6 The lune 

Having extracted everything we can from conformal invariances, we now turn to 
specific manifolds. Evaluation of the zeta function on lunes will turn out to give 
us some more information. Lunes are particularly helpful since the angle between 
adjacent boundary parts is arbitrary. To define a d-lune, we start with a 2-lune, 
which is the region < <fi < 9 on a 2-sphere of unit radius, being the azimuthal 
angle. Higher- dimensional lunes are then defined inductively by ds^_ lunc = dx 2 + 
sin 2 X^(d-i)-i un e> < X < 7T- 
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In our calculations we consider only the case 9 = n/q, where q is an integer, 
and analytically continue our results to all 9. The calculational procedure can be 
extended to 9 = pir/q, with a great deal of added complexity We are assuming that 
the spectral coefficients are smooth functions of 9 - this is the case for C[ d) (0). 



Some work has been done by Cook [14] on the case where 9 /it is irrational. The 
results agree with analytic continuation. 

We will consider 2 < d < 5. Since «Sy ( |14"D is quadratic in d on a spherical 
domain, we can derive at most three independent constraints, and stopping at d = 4 
turns out to be sufficient, d = 5 provides a check of the equations. The procedure is 
to use the eigenvalues and degeneracies to calculate the zeta function and derive the 
heat-kernel coefficients via (ffl). The zeta function on orbifolded spheres, of which 



the lune is an example, has been calculated ]nj, and we do not go into detail here. 
For a general coupling £, we find 



d 2) 



d 4 > 



d 5) 



7T 



720g 
7T 2 a 2 
~3q~ 

Tl 2 



(8g 4 + 20q 2 



7 



+ 



7r 2 a 4 



360g 

vr 3 , 

(11 

45g v 



51 - 60g 2 



2q 

-8g 4 



7ra 



6f/ 



2g 2 



l + 



7ra 



+ 



vr 2 a 2 



3q 



7i 3 a 2 



3q 



2q 2 



2^,4 



3q 



(f-2) 



+ 



7r 3 a 4 



4q 



(35) 
(36) 
(37) 
(38) 



where 



a 



d{d- 1) 



4rf 



In using this data to obtain information about <Sjj, we need to separate off the 
codimension-2 contribution to C% . On the lune all extrinsic curvatures vanish, as 
does V n -R, so the boundary contribution is zero. Since the volume part is propor- 
tional to 1/q, we have 



(<?) 



dp 



i 

—i 

Q 



-a 



The boundary intersection X on a lune is a (d — 2) -sphere of unit radius (2 points, 
each of content 1, in the case d = 2). The boundaries are (d — l)-hemispheres of 
unit radius. The only nonvanishing geometrical quantities are 



R = d(d-1) Ri = Rj = {d - l)(d - 2) R = (d - 2){d - 3) 
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R^n^n^ and R^rijUj vanish since the lune is an Einstein space R^ v = (d — l)g^ v 
and the two normals are orthogonal. 



Equating equation (|I4 ) with our expressions for C 2 
the constraints 



(d) 



on the rf-lune, we obtain 



a i 

3ai + 2a 2 
6a! + 6a 2 + 03 



7T 
12 



(6^-i; 





7T 



7T 



360 



7T" 
03 



e_ 

7T 



7T 




6 V S 



1} 



7r\ 7T 

.tt ~ ?j + 360 





7T 



7T 
03 



from 2, 3 and 4 dimensions respectively. In each case we can conveniently replace 
a\ by Si and remove the (6£ — 1) term. We now have a\, a 2 and for all 0: 

r 3 



a 1 



«2 



a 3 



7T 

360 { 

7T 

~240 
180 1 



7T 
03 





7T 



/ 7T" y 

1^03 ~ n / 
03 TT/ 



(39) 
(40) 
(41) 



^From 5 dimensions we get 

10ai + 12a 2 + 3a 3 



TT 



180 



7T 
03 



TV 



as a check of (p9|), (f|y) and (f41|) . 

The conformal variation equations ( p6|) and ([27|) give 3ai + 2a 2 = 0, in agreement 
with the above. In addition, the lune results agree with equation fl32|) . This is 
encouraging, although we have less information than we hoped for, i.e. lunes have 
only given us 1 new constraint. 

We now have all the a coefficients for general 0, since from (p7|), (p8[) 

1 /TT 3 0\ 

a 4 = ™ + 7^ h?-- cot0 ( 42 ) 



7 A further constraint 

We can in fact make the situation a bit simpler by noticing that h\ is conformally 
invariant for all d, so we can consistently set this quantity to zero. This restriction 
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makes sense since k % 2 vanishes on most spaces we are liable to deal with, e.g. the 
hemiball, lune and cylinder. k\ measures the rate at which rii rotates in the direction 
of Hi as we move around on Zy - an example of a boundary where h\ does not vanish 
is a rectangle twisted so that each edge forms a helix. 

If we set k\ = k\ = 0, there are now only 13 degrees of freedom, and we can 
remove any one of the terms with b or c coefficients from fll4]). We choose to get rid 
of Kikj and KjKi, so we set 65 to zero. 

In the case where this quantity did not vanish, due to its conformal invariance 
in all dimensions we would not be able to find the extra degree of freedom in the 
coefficients except by calculating the zeta function on a manifold where the measure 
of "twist" we have described is non-zero. As far as we can see, this would present a 
very difficult task and for the moment it is necessary to make this constraint if we 
wish to complete the calculation. 

8 Right-angled edges and product manifolds 

We can go further if we limit ourselves to the case 9 = n/2. Again, many of the 
manifolds we come across have this property. We also maintain the constraint we 
have made above, that k\ = 0. For 9 = 7r/2, this condition takes the simple form 

(«i — Kj) 2 = tiK 2 — tr«| (43) 

Our 11 equations, now for only 13 coefficients, become 

_ 7T 7T 7T 1 
Oi — Go = Gq — CLa — 

48 32 24 90 

36i + 6 2 = ^ 263 + 64 = -^ 6i + 56 2 -86 3 = 7 
lo 16 4 

2 , 1 

2ci + 3c 2 = -— — 3ci + 2c 3 + c 4 = d\ = -—— 
135 135 
2 

9ci + 5c 2 + 4c 3 = — (44) 
135 

We can derive the last two constraints we need by considering a product manifold 
M = Mi x M 2 , with an operator A = - (V? + V^) + £(#1 + R 2 ). Mi and M 2 
have smooth boundaries. M has a codimension-2 submanifold dMi x dM 2 with 

9 = tt/2, and obeys flEJ). 

^From (|I|), it is easily shown that, for scalar fields with Dirichlet conditions 

b 2 — O b 2 + Ly l/ 2 Ly 3/2 1 b l + b 3/2 l/2 + b 2 b l 40 J 
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The codimension-2 part of the left hand side is, in terms of geometrical quantities 
on Mi and M 2 



8MixdM2 



d d -' l x\jh x h 2 \ [a x + a 2 )(Ri + R 2 ) + (o 2 + 03) {%. + R 
+ (b 1 + 63) (k\ + nfj + (b 2 + 64) (tm 2 + tm 2 ) + (2ci + c 2 + c 3 )kiK 2 
The expressions for Cjf^ for fc = 0, 1/2, 1, 3/2 are well-known [11, [16 
d d Xy/g 

2 JdM 



a 



(d) 



/nr(«0 
°l/2 



\6 / JM 3 JOM 



°3/2 - 1Q2 



d^xVh \QtiK 2 - 3k 2 + 12(8f - l)i? - 4.R 



(46) 
(47) 
(48) 
(49) 



The codimension-2 part of the right-hand side of fl45|) then comes from the middle 
three terms, and is 



d d - 2 x\lh x h 



dMixdM.2 



1 7T 
— K-i Ko 

9 384 



6 (tr«i + 



-3 (k? + ^) + 12(8f - l){Ri + R 2 ) - 4 (fli + A 
Equating the two sides, we find 



ai + a 2 = — 



71 

96 

7T 



a 2 + a 3 



7T 

96 



61 + 63 = 128 &2 + &4 



7T 

64 



2ci + c 2 + c 3 



9 



(50) 



The first two of these agree with (f44). The second two give one new constraint, 
completing the information for the fe's, and one in agreement with the rest. The last 
completes the information for the c's. Our final result is 



24(1 - 4£)i2 - 12 ( + Rj) + 16R 



Y [ d d - 2 xJj{ — 
faJzv V7 \384 

+24 (trre 2 + tr/? 2 ) + 3 (X 2 + nfj - 30 (tr£ 2 + tr£ 2 

3i?^ (nfn* + n/jnf) - 344 («<«< + Kj Kj) + 228«i« 3 - 



+490^% + 52trKi% - 2 + V~ k 



(51) 
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We note that this, together with equation (]?]), can be extended to the non-zero mass 
case, where the operator becomes 

A = -V 2 + £R - m 2 

simply by replacing £R by — m 2 . 



9 The smeared coefficient and cocycle function 



For completeness, we give the smeared coefficient with 6 = tt/2, where ( f43|) applies. 
The expression for the smooth boundary case is well-known. Using (0), (|5T|) and 
(||), we find 



cf[/] 



M 



d d x^g 



+y [ d^xVh 

i JdMi 



T/+ K«-S) W » / + T5 KV2/ 



— K 2 V n f + — tiK 2 V n f - — V n V 2 / 
210 42 12 . 



+ W d d ~ 2 x^ Ujf + ^-U ( Kl V ni + KjV nj ) f 

(ij) JX v I 04 L 

-9 (^V n , + %V ni ) / + 4 (« + n^nfj 

119 (KiV nj + KjVrn) f - 146 (/? 4 V ni + %V n .) / 



270 



(52) 



where U and T are the volume and boundary parts of the standard smooth expres- 
sion for Cf } [l], ©. 

We can use this to calculate the change in the effective action, or cocycle function, 
in 4 dimensions for a conformal field theory = 1/6). It can be shown using zeta 
function techniques and W = — |C'(0) that under a small conformal variation, 



5W 



1 



The difference in effective action for two metrics g and ~g = e 2uJ g is then given by 

1 n 



W 



W[g\ 



C 



id) 



(4vr) d / 2 Ju=o d ' 2 



-2uw 



g; udu 



(53) 
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Again the volume and boundary parts are known [I7|, the result being 

1 



(W[g]-W[g] 



v,b 2880vr 2 J M 



d x^fg 



^RTupUv + Au^u^u + 2{u> x uj fM ) 2 + 3(V 
1 



2 u) 2 



5760tt 5 



Y { d 3 xVh 

i JdMi 



360Tcj - 2uV n R 



4 16 
--K(V n uj) 2 - —(V n uj) 3 - 2AkV 2 uj - Ak^oo^ 

-20/^cj m + 30V„(V 2 cj + w^Wp 



The calculation of the edge part is fairly easy, and we find 

1 



(54) 



(W[g]-W[g] 



71 

64 



4 [KiVm + KjVnAu 



-9 (^V„. + %V ni ) w + 4 (« + n£n?) 



+3 (V n ^) 2 + (V nj cj) 



- -ujV 2 uj 
3 



119 UiV nj + KjVm) OJ - 146 (a^V^ + %V r 
-2<t#v - 5 + «^) ^ - 58V ni o;V n .a 



a; 



>(55) 



This formula can be used to find the effective action on a manifold if the value on 
a manifold to which it is conformally related is known. Previously, in 4 dimensions, 
we could only consider spaces with smooth boundaries - we can now do calculations 
where an edge is present. 

For example, a hemicap of a 4-sphere of unit radius (~g) is related to the hemiball 
(g) by the conformal transformation 



I /IV 



3 2a ^, u = In f J 



r being the radial coordinate on the hemiball, where the radius a of the hemicap 
and the colatitude of the hemicap boundary are related by 



a = tan — 
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If a = 1, then the hemicap is a quarter-sphere. We find 

W^hemiball( a ) W^l/4-sphere - 



19 , 1 6751 

In a In 2 

1440 360 241920 



W hcmicBp (9) - W 1/4 . 



sphere 



± ln(1 + cose) + _L 



19 In tan 9/2 



Pi 1399 

+- cos 6 (139 + 63 cos 9) cos 3 9 

8 V ; 168 



(56) 



Expression for the quartersphere effective action have been derived in J|, ||. 

10 Spin-1/2 

A similar calculation should be possible if A is the squared Dirac operator 

A = (*7 M V M ) 2 = -g^V,V u + ~R 
as long as the boundary conditions are local. A suitable set are the mixed conditions 

1 



m 19 



0. 



dM V ' n 2 ' 



K tjj 



dM 







for the projection operators P± — |(1 ± ry 5 7 M n M ). 



In this case the volume and boundary parts of C% are known p0|, pl|, p2|. <Sj,- has 



the form (II), and the equations we derive from conformal invariance in 4 dimensions 
are identical to (^) to (^8|), with the exception that 

2d + 3 



Vd 



180 



-V 



where T> is the dimension of the spinor space, and a\ = a\ since there is no variable 
coupling to take into account. 

^From direct calculation on lunes we find 





nV 




720g 




= 




( 




360g V 


cf 






360g V 



(V + 10g 2 + 7 



(57) 
(58) 
(59) 

(60) 
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so that (from the first 3 of these) 



T\V ( 7T 3 7T 6\ 

« = "2880 (V + 10 ? " 17 ;j < 61 > 

llV ( 7T 3 71 e\ 

^ = 192o(V + 10 ?- 17 ^j (62) 

«3 = Z-Ta ( 63 ) 



1440 V 7T # 3 



^From 5 dimensions we get 



TlV ( 7T 3 7T 0\ 

10 ai + 12 a2 + 3a 3 = — (7^3 +40^-47- J 

in agreement with the other relations. As for the scalar case, we have 3ai + 2a2 = 0, 
in agreement with the conformal equations. Additionally, 

IIP TlV ( 7T 2 7T 9\ 

a 4 = 7— + 10- - 17- cot 64 

360 960 V 2 6 tt) K ' 

Unfortunately, this is as far as we can go with the techniques used in this paper. 
Conformal invariance of the smeared coefficient in 6 dimensions does not apply - 
this is to do with the fact that the squared Dirac operator is not in fact conformally 
covariant; it is only a power of a conformally covariant operator. Additionally, on 
the product of two manifolds with boundaries, the heat-kernel does not turn out 
to be a simple product, i.e. equation (|45|) does not apply. This is a result of the 
boundary conditions. 

11 Robin boundary conditions 

A conformally invariant boundary condition, alternative to the Dirichlet case, is the 
Robin condition 

(V n <f>-i/><f>)\ gM = Q (65) 



for some boundary function ip, where under a conformal transformation we set 



so that (|5"5j) is conformally invariant. 



^-i(rf-2)V n cu 



(66) 
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A general expression for the integrand of the codimension-2 contribution is 



Sij = S$ + Ml + K j^i) + u 2 + Kji>i) 

+U 3 (i% + ijjfj + Vi (ftitpj + Kjlpi) + V 2 (Kilpi + Kjlpj) 

+v 3 ipiipj + wi (V-ifji + V-ipj) (67) 



where the Neumann (ipi = ipj = 0) expression Sjj has the same form as the Dirichlet 
expression fll4|), although in general the a, b, c and d constants will have different 
values. 

Going through the procedure we have detailed for the Dirichlet case, this time 
sticking with 6 = ir/2, yields 

_ 7T 7T 7T 1 
ai = do = a 3 = Ola = 

48 32 24 90 

, 7T 7T 57T llX 
Oi = bo = — 63 = Oa = 

1 16 8 128 64 

71 7T 7T 

ui = -- m 2 = - u 3 = - v 3 = 4 (68) 
with 7 equations for the remaining 8 constants: 

2ci + 3c 2 + di + vi = -— (69) 
3d + 2c 3 + Ci + v 2 = (70) 

3di + Wl = -— (71) 
45 

3d + 2v 2 + v 3 + Wl = (72) 

9 11 

9ci + 5c 2 + 4c 3 - 3di + -vi + Av 2 + v 3 - -w 1 = — (73) 

2ci + c 2 + c 3 = - (74) 
y 

2 

d+^2 = ~- (75) 

with checks. 

To acquire the final piece of information, we consider the specific case of a square 
with Neumann conditions on three sides and a non-zero boundary function on the 
fourth. While this analysis would be difficult for general -0, it is possible to calculate 
the eigenvalues and zeta function for a small perturbation ijj — ef, e <C 1. Our results 
will be correct to first order in e. This is sufficient since we will be able to obtain 
the coefficient of the V-^.tbi term in C 2 . 

ill < 
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It can be shown that under such a perturbation, in general the eigenvalues be- 
come 

\ = X N + er]j n=\ d d - l xVh\(f)N\ 2 f (76) 

JdM 

where the A at and <pN are the eigenvalues and normalized eigenfunctions for e = 
(Neumann conditions). 

In the case of a square of side 7r, with ip = ef(y) on the x = boundary, 

4 Fa f 2 1 
Vmn = — } / dy f cos ny, rj 0n = -r] mn 

2 1 

^mo = -o / /, ?7oo = -zVmO) m,n>0 
7l z Jo 2 

where the m, n label the eigenfunctions ~ cos mx cos ray with the eigenvalues 
m 2 + n 2 . 

We will find it convenient to write 

B /l\ 

r/ mn = A + — + O ( -r ) 
where, from integration by parts, 

A = - 2 Tdyf, B = -L tffr) - f (0)] 
Then to first order in e, 

COO = Us)~-se £ E ^+^ + °b) (- 2 + 



2 m "oo"lL W 

-2se?7ooCi?(2s + 2) + (nooe)" 



n 2 



The evaluation of the sum is an old, well-known procedure [23] used in calculating 



the zeta function on cylinders - see for example [|], The standard result is 
£ 1 9(n) (m' + -')" = r( " ' V f £ 

O HZ 00 oo /*oo 

'EE/ ^t u - 3/2 y(n)e-" 2 'e- m27r2 / t (77) 
~, ~I Jo 



IW m=ln=1 Jo 



Since we are interested in C*2 = — 47r£(— 1) we wish to evaluate the above expres- 
sion for u = 0. If g(n) is a finite sum of powers of 1/n 2 , the second term is zero, 
since the sum over m and n converges and lim n ^ r( M ) — 1/u. The first picks out 
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the 1/n 2 term from g(n) since lim u ^o Cr(1 + 2m) = l/2w, and (j? is finite for all other 
arguments. We are left with 

C(-l) = Cat(-I) - l -eBit + 2er ?00 Cfl(0) + er/ o 

=► c 2 = cf + e [/' oo - f (o)] 

since Cr(0) = — 1/2- Comparing with (|67|) immediately gives us the required number 

wi = -1 (78) 

We note that there may possibly be functions f(y) for which this analysis does not 
apply since the sum in (|77|) diverges. However it is of course sufficient to consider 
only one special case - certainly our arguments hold for eg. polynomial f(y), where 
the expansion of g(n) terminates. 

Other heat-kernel coefficients may be evaluated in this way. For example, setting 
s = 0, it is the constant term in g(n) which contributes, and 

C(0) = Ov(0)-^4+l 

In the Neumann case there is a zero mode which is not included in the zeta function, 



so that Cf = 4tt(Ctv(0) + 1), C x = 4vrC(0). Therefore 



C x = Cf-2 / d^xVhip 

JdM 



(79) 



a well-known result. 
Finally, we find 



a 



(d) 



a 



(d)N 



- 5 (mipj + Kjtpi) + 3 (kiipi + Kjtpj) 
+12^ - 3 (V- ^ + V-Vi 



+2 hi + 



10) 



( j(d)N 



£ / d d ~ 2 x^ [24(1 - 4£)J2 - 12 + £,) + 16J2 

+24 (kI + k?) + 48 (tm 2 + tm 2 ) + -15 (a? 2 + £ 2 ) - 42 (tr/? 2 + tm 2 
(n£n? + n^nfj - 434 + + 408/^ 
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+490£;% + 52trK i K j + 88 (V~ + V- k. 
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12 Conclusions 



We have not been able to obtain the co dimension- 2 integrand, Sij, in full for general 
9, except when all extrinsic curvatures vanish. Relaxing this restriction seems to 
be the most desirable extension of our work, and would probably involve a complex 
problem in specific-case evaluation. 

Another interesting extension would be to find for general mixed boundary 
conditions in terms of the P± operators, thus completing the calculations started in 
section 10. As we have noted there, with mixed boundary conditions the heat kernel 
on a product manifold is not in general a product. Work on mixed conditions has 
been done in for the smooth boundary case. 

A problem related to ours is that of conical singularities. These exist on, for 
example, a lune with periodic boundary conditions, so that there is a codimension-2 
submanifold, but not one of co dimension- 1. Fursaev PSfl has evaluated the heat- 



kernel coefficients in this situation, and has derived expressions similar to the in- 
trinsic part of ours. 

The perturbative approach implemented in section 11 may be useful elsewhere 
in the ongoing work of determining heat-kernel coefficients. As an alternative to 
perturbing the boundary function, a perturbation in the coordinates of the boundary, 
and hence in k, yields a simple expression for the first-order effect on the eigenvalues, 
possibly enabling terms in the C*. to first-order in k to be calculated. 

Acknowledgement: The authors would like to thank Klaus Kirsten for some 
very helpful discussions. 
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